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Abstract Why is there no dust in elliptical galaxies? Firstly, the universe is rational and the ratio- 
nality is the proportional density distribution in any independent galaxy. Proportion means that the 
density distribution of stars is orderly. For example, there are four giants standing in array. Their 
heights are A, B, C, D, respectively, and A, B stand in the first row from left to right, C, D in the 



■ second row from left to right. Proportion requires that A divided by C be equal to B divided by D. Ra- 
tional structures are smooth 2- (or 3-) dimensional structures, which are either circularly (spherically) 
symmetric with respect to the center point, or bilaterally symmetric. Secondly, spiral galaxy arms are 

■ linearly shaped, which are neither circularly symmetric nor bilaterally symmetric. Therefore, arms are 
not rational structures. Arm patterns exist only in spiral galaxies and they are weak compared with 

' the main disk structure. Therefore, the presence of arm structure is the disturbance to the rational 

^ ' structure. Because arm patterns exist only in spiral galaxies and only spiral galaxies present dust, 

^ ' I come to the critical answer to the above question: any disturbance to rational structure produces 

cosmic dust. From now on we call the disturbance to rational structure the cosmic impulse. Therefore, 
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the universe is originated not only rationally but also impulsively. 
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■ 1 The Origin of Cosmic Dust 

y3 , 1.1 The Power which Governs Everything: Gravitation 

Scientists have fully proved that there exist only four forces among particles: electromagnetic, weak 
nuclear, strong nuclear, and gravitational. The nuclear forces are short-ranged while the electromag- 
netic force is long-ranged. Each of the three has two contradictory aspects of attracting and repelling 
and, generally, has no net effect in the macroscopic world due to offsetting effect. Gravitational force, 
however, has no contradictory aspects. Gravity has the only effect of attraction and, therefore, can 
not offset itself. Because of this, the true origin of natural structure is the gravitational force. 

The origin of natural structure can not be other forces. Modern science has fully proved that inde- 
pendent system of microscopic particles combined by electromagnetic force or nuclear forces inevitably 
moves towards chaotic state rather than orderly one. This is the principle of entropy increase, which 
is well known for scientists. Therefore, if there were no gravitational force then the whole universe 
would be simply uniform gas without structure. However, there exist in the macroscopic world such 
orderly structures large as galaxies and small as stars, planets, plants, animals, and even human be- 
ings. Therefore, varied kinds of macro-world structures result from the struggling of the gravitational 
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force against the electromagnetic and nuclear forces. 

Unfortunately, gravity is very and very weak. For example, it is 0.0000 • • • 00001 (where 40 zeros 
are after the decimal point) times weaker than the electricity. Only the Earth, Moon, Sun and so 
on present gravity. There is no slightest gravity between cars or human bodies. Therefore, human 
beings suffer insurmountable difficulty to experimentally study gravity. Because human beings are 
insignificant, it is impossible to do physical experiments on such macroscopic systems as the solar 
system or galaxies. However, we can use man-made telescopes to take images of large-scale material 
systems such as galaxies. We can analyze the images. 

1.2 The Origin of Galcixies: Rational Gravity 

The usually familiar gravity refers to the force which exerts between Earth and Moon or between Sun 
and Earth. These examples of gravity are the interaction between two bodies. As for the behavior 
of gravity exerted on many bodies, the solar system can not be the example for us to study such 
behavior. However, each galaxy is composed of billions of stars, which demonstrates the gravitational 
interaction among many bodies. Galaxy images show that each galaxy has a center. Star density at 
the center is the highest. From the center outward, the density is smaller and smaller and presents a 
regular pattern, known as galaxy structure. The principle behind the formation of galaxy structures 
is the demonstration of gravitational interaction in many-body systems. 

Then, what is the behavior of gravitational interaction in many-body systems? we pioneered the 
study on galaxy structures (He (2003, 2004, and 2006), He & Yang (2006)) and the study shows that 
stars in any galaxy are controlled by a very simple orderly force involving many-bodies: proportion. 
Because solar system is just a point at the Milky Way galaxy, the proportion force reduces to Newtonian 
gravity between two bodies. 

Proportion means that the distribution of matters in the universe is orderly. For example, there 
are four giants standing in array. Their heights are A, B, C, D, respectively, and A, B stand in the first 
row from left to right, C, D in the second row from left to right. According to the common view, the 
four giants can have any heights and can stand at any position. That is why current scientific theories 
are incompetent and can not explain the origin of natural structures. They can not provide any basic 
principle to resolve the motion of the most simple gravitational systems (such as interactional free 
three-bodies). However, the universe is orderly. The orderly force at the largest scales requires that 
the distribution of heights is in proportion. In other words, A divided by B is equal to C divided by D. 
This means that A divided by C is equal to B divided by D. If there are nine giants standing in array, 
then the ratios of heights from neighboring two rows are constant (proportion rows). Similarly, the 
ratios of heights from neighboring two lines are constant (proportion lines). In this way are galaxies 
created. 

The above-said rows and lines are all straight (proportion lines). But each galaxy is a regional 
distribution of matters in the universe. Therefore, the proportion lines of each galaxy are curved but 
the rows and columns still cross at each other vertically and they form the net of orthogonal curves. 

In general, a distribution of similar bodies is called the rational structure if its density varies 
proportionally along some particular net of orthogonal curves. Therefore, independent galaxies are all 
rational structures. The force which leads to the rational structure is called the rational force, i.e., the 
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proportion force which is the demonstration of gravity at large-scale and many-body system. That is, 
the universe is rational. This is an important discovery in science. 



2 Two Examples of Rational Structure 
2.1 Logarithmic Density of Galctxy Structure 

Now we start the scientific and mathematical investigation into galaxy structures. A galaxy is a 
distribution of stars. But we can not see individual stars on a galaxy image. A galaxy image is the 
distribution of star densities. Therefore, we use a mathematical function to describe a distribution of 
densities. Because spiral galaxies are planar, we use a function of two variables, x, y, to describe the 
stellar distribution of a spiral galaxy: 

p{x,y) (1) 

where x, y is the rectangular Cartesian coordinates on the spiral galaxy plane. The coordinate origin 
is the galaxy center. Therefore, p{0, 0) is the stellar density at the galaxy center. We want to study 
the ratio of the density p2 to the density pi at two positions 2 and 1 respectively: 

P2/P1. (2) 

In fact, the logarithm of the ratio divided by the distance s between the two positions is approximately 
the directional derivative of the logarithmic density {f{x,y) = In p{x,y)) along the direction of the 
two positions: 

df 

(ln(p2/pi))/s = (lnp2 - lnpi)/s ^ (3) 

There is no systematic mathematical theory on ratios. Therefore, we from now on, study the loga- 
rithmic function f{x,y) instead of the density function p{x,y): 

f{x,y) = lrLp{x,y). (4) 



2.2 Description of a Net of Orthogonal Curves 

The following equation 

x = x{X,fi) 

y = y{\N 

tells us how to describe a net of curves by employing mathematics. Given two functions, x{X,iJ,), 
y{X,IJ,), you have the transformation between the curvilinear coordinates (A,p) and the rectangular 
Cartesian coordinates {x,y), i.e., the formula (5). It describe a net of curves. Letting the second 
parameter p be a constant, you have a curve (called a row curve, i.e., the proportion row defined in 
the following). That is, the formula (5) is a curve with its parameter being A. For the different values 
of the constant p, you have a set of "parallel" rows. On the other hand, you let the first parameter A 
be a constant then you have another curve (called a column curve). That is, the formula (5) is a curve 
with its parameter being p. For the different values of the constant A, you have a set of "parallel" 
columns. 
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However, The row curves and the column curves are not necessarily orthogonal to each other. The 
following equation is the necessary and sufficient condition for the net of curves to be orthogonal: 

dxdj^ dy_dy_^ 

dXdfi dXdfi~ ' ^ ^ 

To study the rational structures described in the following, we need more knowledge of the de- 
scription of row and column curves. The arc length of the row curve is s whose differential is the 
following: 



ds = ^x'^ + y'^dX = PdX. (7) 
where P is the arc derivative of the row curve: 



P = A = ^|x'l + y'l (8) 

The arc length of the column curve is t whose differential is the following: 

dt = yjx'^ + y'^dii = Qd^i. (9) 
where Q is the arc derivative of the column curve: 

Q = ^; = ^^i^^y^F■ (10) 

2.3 The Condition of Rational Structure 

The formulas (5) and (6) are the general description of a net of orthogonal curves, and the formula (1) 
is the general description of a distribution of densities (a structure). This paper talks about rational 
structure. A distribution of densities is called the rational structure if its density varies proportionally 
along some particular net of orthogonal curves. That is, you walk along a curve from the net and the 
ratio of the density on your left side to the immediate density on your right side is constant along the 
curve. However, the constant ratio of this curve is generally different from the constant ratios of the 
other curves. 

We have shown that a logarithmic ratio of two densities divided by the distance between the two 
positions is approximately the directional derivative of the logarithmic density along the direction of 
the distance. Therefore, we always study the logarithmic function f{x,y) (see the formula (4)). If we 
know the two partial derivatives df/dx and df/dy then the structure f{x,y) is found. The partial 
derivatives of f{x, y) are the directional derivatives along the straight directions of the rectangular 
Cartesian coordinate lines. However, we are interested in the net of orthogonal curves and what we 
look for is the directional derivatives along the tangent directions of the curvilinear rows and columns. 
These are denoted u{X,fi) and t;(A,/i) respectively. 

The condition of rational structure is that u depends only on A while v depends only on //: 

V = V{fJ,). 

What a simple condition for the solution of rational structure! 

Now we prove the condition. Assume you walk along a row curve. The logarithmic ratio of the 
density on your left side to the immediate density on your right side is approximately the directional 
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derivative of f{x,y) along the column direction. That is, the logarithmic ratio is approximately the 
directional derivative t;(A,/i). Because v{\, fi) is constant along the row curve (rational), v(X,jj,) is 
independent of A: v = v{iJi). Similarly, we can prove that u{X,iJ,) = u{X). 

2.4 The Equation of Rational Structure 

It is known that, given an arbitrary function, we can have its two partial derivatives. However, given 
two functions, we may not find the third function whose partial derivatives are the given functions. 
For the given two functions to be some derivatives, a condition must be satisfied. The condition is 
the Green's theorem. In the case of derivatives along orthogonal curves, the Green's theorem is the 
following: 

|-(n(A,/x)P)-|^(KA,M)<3)=0 (12) 

In the case of rational structure, directional derivatives arc the functions of the single variables, A 
and n respectively (see (11)). Therefore, the Green's theorem turns out to be the following which is 
called the rational structure equation: 

u{X)P'^-viMx = (13) 

This equation determines rational structure. To find a rational structure, generally we are first of 
all given a net of orthogonal curves. Accordingly the arc derivatives of both the rows and columns, 
P(A, /x) and Q(A, /x), are known. Therefore, the remaining functions , u{X) and v{fi), are the only 
unknowns. Note that the unknowns are the functions of the single variables, A and fi respectively. 
Because the rational structure equation involves no derivative of the unknowns, the equation is not a 
differential equation at all. It is an algebraic equation and what we need to do is to add two factors 
(the two unknowns) to the derivatives of P and Q so that the rational structure equation holds: 
factorization! 

How simple is the universe constructed! 



2.5 The First Example: Spiral Galaxy Disk Structure 

Images of spiral galaxies taken with infrared light show that each spiral galaxy is mainly a disk (the 
circularly symmetric disk with respect to the galaxy center, i.e., the disk center) and the disk light 
density decreases exponentially outwards along the radial direction from the center. There are other 
minor or weak structures in spiral galaxies. Spiral galaxies gain their name by the fact that they 
present more or less spiral structures, known as arms. 

The first example of rational structure is the spiral galaxy disk and is determined by the following 
net of equiangular spirals (or called logarithmic spirals) (see Fig. 1): 

X = e''i^+'^2^ cos((i3A + d4,fi), 

y = sin(d3A + d^n) ^ ' 

where di(> 0),d2(> 0),o?3(< 0),o?4(> 0) are real constants and we choose 

ds = -did2/d4 (15) 
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so that the curves are orthogonal. The polar angle and polar distance of the point {x,y) are easily 
seen, 

r — f>diX+d2lJ- 

(16) 

9 = d^X + d4ii. 

The coordinate lines are spiral-shaped and shown in Fig.,1. 
The corresponding arc-length derivatives P, Q are 



P(A, ^) = s', = {di/d,)^dl + dle'^l^+'^2^ 

Q{\y.) = t'^ = ^di + dy^^+''-'^. 

The rational structure equation (13) does help factor out the required directional derivatives for our 
spiral galaxy disks, 

u{\) = d^di, ^^g^ 
t;(/x) = 

where d^ is another constant. 

So far, we did not specify the variance domain S on (A, /x) coordinate plane on which the coordinate 
system (14) is defined. There are many such domains which are proved in the Appendix of the paper. 
Here we simply present one kind of such domain. We define a constant. 

The domain is the following: 

Sai : Ai < A < A2(= Ai + Aa), -oo < /x < +oo. (20) 

where Ai is arbitrary constant and the length of the interval (Ai, A2) is 

A2 = Ai + Aa. (21) 

In the Appendix, we also prove that the spiral curves are equiangular spirals. 

The spiral disk density along all orthogonal coordinate lines can be found by performing path 
integrations of the following formulas 

df = uds = u{X)PdX, ^22^ 
df = vdt = v{jjL)Qd^. 

along the row curves, /x = constant, and the column curves, A = constant, respectively. Without loss 
of generality, we assume the coordinates A and /U are defined on the domain (20). Then starting at 
the galaxy center, we perform the path integration of df = v{iJ,)Qdfi over [— oo,;[i] by taking A to be 
the constant Ai to get /(Ai,/i). Then we perform the path integration of df = u{X)PdX over [Ai,A] 
by taking fj, to be an arbitrary constant to get /(A, /n). Finally, we have the logarithmic function 
f{x,y) implied by the spiral-shaped coordinate system. The density distribution p(x,y) represents 
spiral galaxy disks (we choose ^5 < 0, because light density p — >■ when r — >■ -|-oo). 



fd = d,^d?, + dy^^+^^'^, 

Pd = doexp(4v^i + d|e«'i^+'^2/.) 
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where do is the light density (star density) at the galaxy center. Note that we use the letter d as well 
as the subscript d for disk parameters and formulas. Similar notations are used for bar parameters 
and formulas. Because fa = d^Q{\,^) and Q{X,fi) is uniquely defined over the whole galaxy disk 
plane ({x,y) plane, see the Appendix A), fd{x,y) is uniquely defined on the same plane. We can see 
the disk light pattern by displaying pii{x,y). Because the polar distance is r = exp(ciiA + ci2yu), the 
galaxy disk light distribution is circularly symmetric, 

= doe^'^ = doe(<^5V^I+rf!)r_ (24) 

Since galaxy light density oc p, we have recovered the known exponential law of spiral galaxy disk (the 
exponential disk). 

The circularly symmetric exponential disk is completely determined by the value of do the 
value of d5'\J(i2 + d\. Therefore, given a spiral galaxy disk, that is, given the two specific values, we 
can find an infinite number of coordinates (14) defined on the disk plane which give the same disk 
structure. From now on, we choose y^dl + ^4 = 1 so that the disk model involves only two variable 
disk parameters (io(> 0) and d^{< 0), 

Pd = doe^^ = d^e"^^'. (25) 

In fact, galaxy structures depend only on the geometric curves, not the choice of coordinate 
parameters. Therefore, the following is the general expression for equiangular spirals: 

X = e'^i/(^)+'^29(/.) cos(d3/(A) + d^g{ii)), 
y = e'^im+d^giu.) sin(d3/(A) + d4g{p)) 

where f{\),g{p) are arbitrary functions. All these expressions give the same equiangular spirals and 
generate the same exponential disks, independent of the choice of coordinate parameters. That is, 
our rational structure which depends only on the geometric proportion curves, is independent of the 
coordinate system which describes the curves (coordinate invariance). 
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2.6 The Second Example: Dual Handle Structure 

Now we study galactic bar model. A bar pattern is composed of two or three dual handle structures 
which are generally aligned with each other (spiral galaxy NGC 1365 is not the case). The dual handle 
structure is determined by the following orthogonal curves of confocal ellipses and hyperbolas: 

X = e'^ cos T, 



y = ^e^" + bl sinr, (27) 
—GO < a < +00, < r < 27r. 

where 6i(> 0) is a constant. The orthogonal coordinate system no longer shares the coordinate lines 
with the polar coordinate system. The coordinate lines are confocal ellipses and hyperbolas (Fig. 2). 
The distance between the two foci is 2bi which measures the distance between the two handles. The 
eccentric anomaly of the ellipses is r. 

Now we look for the logarithmic density (/(cr, r)) of the dual handle structure determined by the 
coordinate system (27). The corresponding arc-length derivatives of the coordinate system are (see 
(8) and (10)), 

P = e-^e2- + 6fcos2r/<^e2- + 6f, 
Q = ^e2<^ + 6fcos2r. 

The rational structure equation 

u,{a)P^ - v,{t)Q', = (29) 
determines the corresponding directional derivatives of the logarithmic density, 



ui,{a) = 626°^-^ e2<^ + h\, vi,{t) = —6261 sinr cost. (30) 

To get the logarithmic density we need to perform path integrations of the similar formulas to (22). 
The result is 

h{cT, r) = (62/3)(e2- + bj cos^ rf. (31) 
The inverse coordinate transformation of the formulas (27) is easily found, 



p{x, y) = e- = ^(r2-6f + 7(72-62)2 + 462x2)72, ^g^) 
COST = xe~^ = x/p{x, y) 

where r'^ = x^ + y^. Finally we find the density of the dual handle structure, 

Mx,y) = {b2/3){pHx,y) + bW/pHx,y))^/^, ^^^^ 
Pb = boexp{fb{x,y)) 

where bo is the dual handle density at the galaxy center. We need to choose 62 < so that fb<0 
and — > when r — > +00. We can see that 60 corresponds to the central dual handle strength and 
61 corresponds to the dual handle length while 62 measures the density slope off the dual handle. If 
we display the dual handle structure as a curved surface in 3- dimensional space then we can see that 
the surface is camelback-like shapes with two humps (i.e., handles). 
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Figure 1: The orthogonal net of logarithmic spirals (14). The angle i, at each position on the spiral, 
between its bending direction and the disk radial direction is constant along the curve. The closed 
curve which consists of two sections of equiangular spirals (thick dotted lines), demonstrates the closure 
condition (see the Appendix A for details). 
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Note that polar angle 9 is not defined at the center r = for the polar coordinates. Similarly the 
eccentric anomaly r is not defined for the coordinates (27) along the dual handle central line of 2b\ 
length which is the coordinate line a = — oo. 

In fact, galaxy structures depend only on the geometric curves, not the choice of coordinate 
parameters. Therefore, the following is the general expression for dual handle structures: 



where /(A),5f(/x) are arbitrary functions. All these expressions give the same dual handle structures, 
independent of the choice of coordinate parameters. That is, our rational structure which depends 
only on the geometric proportion curves, is independent of the coordinate system which describes the 
curves (coordinate invariance). 




(34) 
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3 The Origin of Spiral Galaxies 



3.1 Proposition!: Rational Structures are at Most Bilaterally Symmetric 

This is a mathematical proposition: any net of orthogonal curves is either circularly symmetric with 
respect to the center point, or bilaterally symmetric. I spent three years from 2002 to 2005 to look 
for a net of orthogonal curves whose shape has odd symmetry. That is, I wanted to find a rational 
structure which resembled a two-arm spiral pattern like the spiral galaxy M51. The three-year study 
indicates that a net of orthogonal curves is generally connected to some complex analytical function. 
As you might know, a complex analytical function always has two parts (real and imaginary) like the 
formula (5), and leads to a net of orthogonal curves. However, Complex analytical functions are very 
special ones which satisfy some strong conditions like Cauchy integral theorem and formula. Prom my 
experience, I do not find any complex analytical function whose graph of the real or imaginary part 
has odd symmetry. Therefore, I have the proposition that rational structures are at most bilaterally 
symmetric. This is left for verification by able people. 

Surprisingly, galaxy structures happen in the same way as indicated in the following. 

3.2 Coincidence!: Galstxy Patterns (except Arms) aire at Most Bilaterally Sym- 
metric 

Amazingly, any component of any galaxy pattern (except the arm pattern) is either circularly sym- 
metric, or bilaterally symmetric. And my academic papers (He, 2005a, 2005b, and 2008) show that, 
except the arm structure, any component of any galaxy structure (such as exponential disks, galactic 
bars, even the whole elliptical galaxy) is a rational structure. That is, any galaxy image (ignoring the 
arm) can be fitted identically to rational structures. 

The following will present more and more cases of coincidences. Therefore, they are not coincidence 
at all. They are the cosmic truth, and my discovery is absolutely important. 

3.3 Coincidence 2: Dust and Irrationality 

It is the observational fact that spiral galaxies are full of dust while elliptical galaxies have no dust. 
Why does this happen? However, people have not found the appropriate answer since galaxies were 
discovered more than 80 years ago. 

Arm structure is neither circularly symmetric with respect to the center point, nor bilaterally 
symmetric. Therefore, arms are not rational structures. Arm pattern tends to be odd symmetry with 
respect to the center. But we can never find such "grand design" arm pattern which presents the 
perfectly odd symmetry. In fact, there are very different types of spiral structures. Some galaxies are 
what we call "Grand Design" spirals, meaning that they have a clearly outlined and well organized 
spiral structure. Other galaxies, like NGC 4414 are called "fiocculent" and have much harder to trace 
arms. Arms may not be called a structure at all as compared with the exponential disks, galactic 
bars, and even the whole elliptical galaxies which demonstrate smooth and rational structures, . 

Arm patterns exist only in spiral galaxies and they are weak compared with the main disk structure 
of spiral galaxies. Therefore, the presence of arm structure is the disturbance to the rational structure. 
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Figure 2: The orthogonal curves of confocal elhpses and hyperbolas. The distance between the two 
foci, F and F', is 2bi which measures the distance between the two handles. 
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Because arm patterns exist only in spiral galaxies and only spiral galaxies present dust, I comes to the 
critical answer to the above question: any disturbance to rational structure produces cosmic dust. 

3.4 The Origin of Cosmic Dust: Impulsive Gravity 

Prom now on we call the disturbance to rational structure the cosmic impulse. Therefore, the universe 
is originated not only rationally but also impulsively. In the case of large-scale structures (that is, 
galaxies), the impulse is demonstrated to be the disturbing waves, i.e., the arm patterns. 

We have seen that the rational structure is the main structure while the disturbance to the main 
structure is always weak. That is, rational force is the main one while the impulse is weak. In the case 
of large-scale structures (that is, galaxies), the disturbing waves try to achieve the minimal disturbance 
and, as a result, they follow the proportion rows or columns of the rational structures. The impulsive 
disturbance reveals the rational design of the universe: proportion. 

Because we live inside a galaxy (Milky Way galaxy) and wc find no other force exists except gravity 
at our neighborhood, we conclude that Newton or Einstein formulation of gravity is a partial result 
of the universal "gravity" that not only presents rational sense but also impulsive one as indicated at 
the large-scale galaxy system. 

3.5 A Simpler Argument that Spiral Galeixy Disk is a Rational Structure 

The logarithmic stellar density of exponential disk is the circularly symmetric distribution of numerical 
values about the center point, which decreases linearly in the radial directions. Therefore, the magni- 
tude of the directional derivatives along radial directions is a global constant. An equiangular spiral 
is the curve which makes a constant angle to the local radial directions. Therefore, the magnitude of 
the directional derivatives along the tangent directions of an equiangular spiral is constant along the 
spiral. Accordingly the magnitude of the directional derivatives along the perpendicular directions to 
the equiangular spiral is also constant along the spiral. That is, you walk along an equiangular spiral 
and the ratio of the density on your left side to the immediate density on your right side is constant 
along the curve. However, this constant is different from the constant in the radial direction. They 
differ by a factor of the cosine of the angle. Finally we have proved that exponential disk is a rational 
structure. 

All curves which are orthogonal to a specific set of equiangular spirals in the clockwise direction are 
themselves equiangular but in the counter-clockwise direction. These two sets of spirals in opposite 
directions consist of the orthogonal net of curves, and the exponential disk is the rational structure 
with respect to these proportion curves. 

3.6 Coincidence 3: Exponential Disk is Correlated with Equiangulair Spiral 

The exponential disk of any spiral galaxy is a rational structure which is circularly symmetric about 
the galaxy center. It is amazing that its proportion curves are the equiangular spirals which are 
precisely the curves represented by normal spiral galaxy arms. This is another coincidence. 
Now we introduce more and more coincidences. 
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3.7 Coincidence 4: The Only Rational Structure which is not Circularly Sym- 
metric is Dual Handle Structure 

I have given the definition of rational structure. However, given an arbitrary net of orthogonal curves, 
we are not always possible to arrange a distribution of stellar density on the net to form a rational 
structure. In fact, there are only a few types of orthogonal curves which correspond to rational 
structures. Rational structures arc usually circularly symmetric about the center points. The only 
rational structure we can find which is not circularly symmetric, is the bilaterally symmetric structure, 
namely, dual-handle structure. 

We have two types of rational structures: exponential disks and dual-handle structures. Adding 
the two structures together leads to the barred pattern as we expected. It is amazing that only 
two kinds of spiral galaxies are observed in the universe. One kind of spirals are the normal spiral 
galaxies while the other kind are the barred spiral ones. What is more surprising is that some barred 
galaxies do show a set of symmetric enhancements at the ends of the stellar bar, called the ansae or 
the "handles" of the bar (see upper- left panel of Fig. 3). This indicates that a bar itself is nothing 
but a set of several pairs of ansae (handles). That is, bars are the superposition of several aligned 
or misaligned dual-handle structures. If the outer dual-handle structure is far more away from the 
galaxy center then it demonstrates the pattern of ansae or "handles" of the bar. 
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3.8 Coincidence 5: There are Barred Spiral Galcixies which Present Two Non- 
parallel Bars 

The main structure of spiral galaxies is the exponential disk. When the dual-handle structure (i.e., 
sub-bar) is near the galaxy center, the superposition of the dual-handles to the bright disk center 
presents a bar shape. This precisely explains the origin of galaxy bars. A galaxy bar is usually 
composed of two or more sets of aligned or misaligned dual-handle structures. Surprisingly, there are 
barred spiral galaxies which present two nonparallel bars (see upper-right panel of Fig. 3). 

3.9 Coincidence 6: Bcir Structure is So WeEik in the Outer Areeis of Spiral GalzLxies 
that it is Ignored 

Compared with the exponential disk, the bar is observationally weak structure. That is, bar structure 
is so weak in the outer areas of spiral galaxies that it is ignored. It is very surprising that the theoretical 
calculation of dual-handle structure shows that it is weak when compared with the disk (see Fig. 4, 5, 
and 6). 

We know that the disk density of spiral galaxies decreases outwards exponentially, which is the 
numerical result obtained over 90 years since the discovery of galaxies in the universe. Spiral galaxy 
disks are thus called exponential disks. We add the dual-handle structure to the exponential disk for 
them to be the model of barred spiral galaxies. If the density of dual-handle structure were comparable 
to or stronger than the exponential disk in the far distances from the galaxy center then our model 
would fail. That would suggest that the main structure of spiral galaxies were not the exponential 
disk, a result inconsistent with astronomical observation. The mathematical result is that the density 
distribution of dual-handle structure decreases outwards cubic-exponentially as shown below. 

For large r, the first formula in (32) reduces to r, and the first formula in (33) approaches (62/3)r^. 
Because 62 < 0, we come to the conclusion that the density distribution of dual-handle structure 
decreases outwards cubic-exponentially. In the distances far from the galaxy center, the dual-handle 
structure is negligible when compared with the disk. Mathematical result is consistent to observation. 
This piece of coincidence alone proves that galaxies are originated from proportion force. 
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Figure 3: Upper-left panel is galaxy NGC 2983 (image credit: Martinez- Valpuesta, Knapen, & Buta 
(2007)). Upper-right is the infrared image of galaxy NGC 1365. Lower-left is the infrared image of 
NGC 1300. Lower-right is NGC 1365 (ultraviolet image, credit: European Southern Observatory). 
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3.10 Coincidence 7: The Arms of Barred GalcLxies Spin ciround the Bar and cire 
No Longer the Equiangular Spirals 

With simple mathematical calculation we know that spiral-shaped proportion curves exist in dual- 
handle structure. However, they are not equiangular because they surround the central line of the 
dual-handles (recall that the spirals in exponential disks are equiangular and surround the center 
point). Two proportion curves which arc oddly symmetric about the center point in dual-handle 
structure make approximately elliptical shape and its long axis must be parallel to the central line 
of the dual-handles. Surprisingly, astronomical observations show that arms of barred spiral galaxies 
do surround the middle lines of their bars, and they are not equiangular, and the two arms make 
approximately elliptical shapes with the long axes being parallel to the bar middle lines (see Fig. 3) . 
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We proceed with more mathematical details. The proportion curves of a dual-handle structure are 
all confocal ellipses and hyperbolas. The two foci are the centers of the two handles. The distance 
between the two foci is the length of the dual-handle structure. 

If the length of the dual-handle structure is zero then the two foci overlap to be the center of 
concentric circles and the above-said proportion curves become the curves of polar coordinates. This 
returns to the case of normal spiral galaxy disk. In fact, the polar curves are also the proportion 
curves of normal spiral galaxies. In other words, all polar curves are the limiting curves of equiangular 
spirals. 

Back to the above-mentioned dual-handle structures. Similarly, they have also open proportion 
spirals which make acute angles to the above-mentioned ellipses and hyperbolas, and spin around the 
central line of the dual-handle structure. However, the spirals are no longer equiangular. 

3.11 Coincidences: Circular and Elliptical Rings 

We have already known that exponential disks have circular proportion curves (one family of polar 
curves). Observationally, some normal spiral galaxies do have closed arms which are circular, called 
rings. Dual-handle structures also have closed proportion curves which are ellipses whose long axes 
must be parallel to the central lines of the dual-handles. Observationally, some barred spiral galaxies 
do have closed rings which are ellipses and the long axes are parallel to the galaxy bars. 
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Figure 5: The OSUBGS H-band images NGC 4665, 4930, 5701 (image credit: Eskridge, et al. (2002)) 
minus our model bars respectively result in the disk and bulge images (bar deletion). 
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3.12 Coincidence 9: Fitting Bar Images with Dual-handle Structures 

In fact, we can directly use dual-handle structures to fit into the bar images of spiral galaxies. Prom 
the barred spiral galaxy images we subtract the superposition of two or three sets of dual-handle 
structures and sec if the resulting images are exponential disks plus the weak arm structures. If it 
works then we have further proved that barred spiral galaxies are the superposition of exponential 
disks with dual-handle structures. The first column of the Figures 4, 5 and 6 present real galaxy 
images and the second column present the fitting bars. The third column presents the result of the 
subtraction of the second column from the first one. The results are very good: the resulting images 
are exponential disks plus the weak arm structures. In fact, I have made a piece of computer software 
which can be used to improve the results. 

3.13 Coincidence 10: Elliptical Galetxies aire Completely Rational Structures 

I have proved that elliptical galaxies are completely rational structures in three-dimensions (He, 2005a 
and 2008). The proportion surfaces of elliptical galaxies are the intersecting nets of orthogonal spheres, 
where disturbing waves are difficult to form and spread. On the other hand, spiral galaxies are two- 
dimensional and their proportion curves are open spirals where disturbance waves are easy to form 
and spread. Astronomical observations do show that arms do not exist in elliptical galaxies. 

The disturbance to rational structure leads to the formation of gas and dust. New families of stars 
and planets are born to these gas and dust. The star-planet families are short-lived. This happens 
only in spiral galaxies. 

3.14 It Seems that Fitting Galaxy Images may Tell the Physical Distances of the 
Corresponding GalEtxies 

Galaxies are very far away from the earth and human is really minute. Therefore, we can not measure 
directly how far away the galaxies arc from Earth. It seems that fitting barred galaxy images with 
dual-handle structures may tell the physical distances of the corresponding galaxies. However, the 
result needs further confirmation with the help of supercomputers. 

3.15 Future Work and Coincidence 11: Galctxy Inner Rings. 

There are lots of future theoretical and applied works to be done. One of them is the search of the 
condition under which the summation of two rational structures is still rational. Let pi and p2 be the 
two rational structures. Their logarithmic densities are /i and /2 respectively. The logarithm of the 
summation is /. It is straightforward to show that 

V/= (piV/l+p2V/2)/(pi+P2). 

That is, the gradient of summation is the summation of weighted gradients. Because galaxy disks 
take their maximum values at around the galaxy center and dual handle structures take the values 
away from the center, our above conclusions are still true. Furthermore, the proportion curves of 
the summation structure explain why there are nuclear rings or arms for some barred galaxies. An 
example is the galaxy NGC 4314. 
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A The Equiangular Coordinate System as a Uniquely Defined Map- 
ping between Coordinate Spaces. 

We did not specify the variance domain S on (A, //) coordinate plane on which the equiangular coor- 
dinate system (14) is defined and maps it onto the whole (x, y) coordinate plane. To find the domain, 
we define two constants which are called periods, 

^-^ didi-d^d-z C35-) 
— d-idi-d3d2 V-^ 

It is straightforward to show that the two periods satisfy the following equations, 

di^x - d2^u = 0, 

-dgAA + d4A^ = 27T. 

There arc many such domains and now wc want to prove the statement of coordinate periodicity 
that any vertical infinite band Sx^ on (A, /i) plane can be mapped onto the whole {x, y) plane of a 
galaxy disk, 

: Ai < A < Ai + Aa(= A2), -00 < /x < +00. (37) 
where Ai is arbitrary constant and the length of the interval (Ai, A2) is the period, 

A2 = Ai + Aa. (38) 

As indicated in Fig. 1, it is equivalent to show that the two different vertical boundary lines, A = Ai 
and A = A2 on (A,//) plane are mapped to a single curvilinear coordinate line on {x,y) plane (see 
Fig. 1). This can be shown by three steps. We choose a closed curve (thick dotted line in Fig. 1) which 
consists of two sections, one from the coordinate lines of first set (/x = constant = Hi,Xi < A < A2) 
and the other from the second set (A = constant = Ai,/xi < n < 112) ■ The closed curve is called 
snail-shaped curve (see Fig. 1). 

For the first step, we show that the arc-length derivative Q{\,fi) is uniquely defined along the 
snail-shaped curve. Starting at the point Ni where the value of Q is (5(Ai,/xi), we have two ways to 
go to the point N2, one being the section F and the other T. The corresponding two values of Q at 
N2 are Q{\i,^2) and Q(A2,^i) respectively. The uniqueness of Q means that (3(Ai,/X2) = <5(A2,A*i). 
This is guaranteed by the first equation in (36). 

For the second step, the values of the tangent angle a(A, //) of the coordinate line A = Ai (see 
Fig. 1) should be uniquely defined along the line within a difference of 27r. Starting at the point Ni 
where the value of a is a{\i,iii), we have two ways to go to the point N2., one being the section F 
and the other T. The corresponding two values of a at N2 are q;(Ai,/X2) and a{\2,Hi) respectively. 
Therefore, 

(q;(Ai,/X2) - a(Ai,/xi)) - (a(A2,/xi) - a(Ai,/xi)) = 2it. (39) 
Now we need calculate the value of a(A, //). Its calculation is straightforward with formulas (14), 

y'a d2 sin 6 + diCos9 

tana = -j = ^ — , . ^ - (40) 

02 cos p — 04 sm u 



23 



Defining 



cos ^0 = 



, sin^o = 



\ldl + dl 



(41) 



then we have 



tana = tan(0 + ^o)- 



(42) 



Finally we find the value of the tangent angle of the coordinate line A = constant 



a = 9 + 9q = + d^jjL + ^o- 



(43) 



Substituting the value into the equation (39), we find that the equation reduces to the second equation 
in (36) which is already proved. Note that, in Fig. 1, a is the tangent angle of the coordinate line 
A = constant at a point {x, y) and 9 is the polar angle of the same point. Their difference is the pitch 
angle of the coordinate line A = constant. Therefore, we find the pitch angle i to be, 



which is constant over the whole plane. We have proved that our coordinate lines (proportion curves) 
are equiangular spirals. 

For the third step, we need prove that the snail-shaped curve is really closed. To see the case, we 
change the coordinate transformation equations (14) into different form, making use of the formulas 



Because Q(A, /x) and a are uniquely defined and both and -^(dl+dl) are constants, x, y are uniquely 
defined. That is, the snail-shaped curve is really closed. 

Therefore, we have shown that the two different vertical boundary lines, A = Ai and A = A2, 
on (A,/x) plane are mapped to a single curvilinear coordinate line on (x,y) plane (see Fig. 1). This is 
equivalent to say that the vertical infinite band on (A, ^) plane (see (36)) is mapped onto the whole 
(,T, y) plane of the galaxy disk. Because Ai is arbitrarily chosen, we have shown that the coordinate 
transformation (14) is periodic in that each interval of A of length (period) A;^ with —00 < jjb < +00 
is mapped onto whole {x, y) plane. The periodicity of can be similarly proved. 



i = a — 9 = 9q 



(44) 



(17) and (43) 




(45) 
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